X=[1 0; 1 1; 1 2; 1 3];
y=[1 1 2 2]';
beta=X\y